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ABSTRACT 

Cosmic shear requires high precision measurement of galaxy shapes in the presence 
of the observational Point Spread Function (PSF) that smears out the image. The 
PSF must therefore be known for each galaxy to a high accuracy. However, for several 
reasons, the PSF is usually wavelength dependent, therefore the differences between 
the spectral energy distribution of the observed objects introduces further complexity. 
In this paper wc investigate the effect of the wavelength-dependence of the PSF, 
focusing on instruments in which the PSF size is dominated by the diffraction-limit 
of the telescope and which use broad-band filters for shape measurement. 

We first calculate biases on cosmological parameter estimation from cosmic shear 
when the stellar PSF is used uncorrected. Using realistic galaxy and star spectral 
energy distributions and populations and a simple three-component circular PSF we 
find that the colour-dependence must be taken into account for the next generation of 
telescopes. We then consider two different methods for removing the effect (i) the use of 
stars of the same colour as the galaxies and (ii) estimation of the galaxy spectral energy 
distribution using multiple colours and using a telescope model for the PSF. We find 
that both of these methods correct the effect to levels below the tolerances required 
for per-cent level measurements of dark energy parameters. Comparison of the two 
methods favours the template-fitting method because its efficiency is less dependent 
on galaxy redshift than the broad-band colour method and takes full advantage of 
deeper photometry. 
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1 INTRODUCTION 

Measurements of the cosmic shear signal are expected to 
play a leading role in furthering our understanding of our 
Universe, in particular the nature of dark matter and dark 
energy or its possible alternatives such as modifications in 
gravity. The gravitational lensing of light from distant galax- 
ies by intervening mass provides a powerful insight into the 
growth of structure an d the expansion history of the universe 
(for recent reviews see iHoekstra fc Jainll2008l : [ Munshi et al.l 
[2008 : RcfrcEicr 2Q0i). 

Several planned future dark energy missions are de- 
signed with weak lensing as a primary science driver. 



including ground-based projects: the Kilo-Degree Survey 
(KIDS) 3 the Panoramic Survey Telescope and Rapid Re- 
sponse System (Pan-STARRS) Q; the Dark Ener gy Survey 
(DES)0; and the Large Synop tic Survey Telescope (LSST)B 
and s pace missions Euclid FT (|Laureiisll2009l : iRefregier et al.l 
|2010|) and the Joint Dark Energy Mission (JDEM) Q. The 
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success of the me thod relies on accurate measurement of 
galaxy shapes (e.g . lHevmans et al ]|2006l : lMasse7 et al ] |2007l : 
iBridle et al.ll2009l ). 

As light from galaxies passes through the atmosphere, 
telescope optics and measurement devices it is convolved 
with a kernel, referred to as the Point Spread Function 
(PSF). The size of the PSF is similar to the size of the galax- 
ies used to measure cosmic shear, making accurate determi- 
nati on of the underlying galaxy shape a significant challenge 
(e.g. lLewisll2009l : IVoigt fc Bridi3l2009l). Precise model ling of 
the PSF is crucial (e.g. IPaulin-Henriksson et al.ll2008l 'l. This 
can in principle be performed using accurate knowledge of 
the telescope optics, however in practice it is usual to make 
use of the point-like nature of stars by modelling the PSF 
based on the shapes of stellar images. The situation is com- 
plicated by (i) variations in the shape of the PSF across 
the detector and (ii) the wavelength-dependence of the PSF 
shape. In this paper we concentrate on the second of these 
two effects, focusing on a telescope which is nearly diffrac- 
tion limited and on a instrument that uses broad-band op- 
tical filters for shape measurement. 

Stars and galaxies have different spectral energy dis- 
tributions (SEDs), both intrinsic and observed, specially if 
we consider that galaxies are observed over a broad range 
of redshifts (See Figure [T|. Therefore the PSF shape mea- 
sured from a stellar image will be different from the true 
PSF applied to the galaxy, leading to a bias on the mea- 
sured underlying galaxy shape. The wavelength-dependence 
of the PSF has a greater impact if observations are per- 
formed using broad-band filters. In this paper we quantify 
the effect of ignoring the wavelength-dependence of the PSF 
in a Euclid-like space mission and show the feasibility of re- 
ducing the effect by using colour information. We determine 
the wavelength-dependence of the PSF shape using a sim- 
ple model for the convolution kernel and applying realistic 
galaxy and stellar spectral energy distributions (SEDs). 

For the reasons we will discuss in the following sections 
this paper will focus on the wavelength-dependence of the 
PSF size. There are several measures of this quantity, such 
as the 50we adopt the full width at half-maximum intensity 
(FWHM) as our size parameter given that this is a more 
familiar quantity to the astronomical community and de- 
scribes well the simple PSF models we use in this paper. 

The paper is organised as follows: in Section [5] we set 
out a formalism for propagating a mis-estimation of the PSF 
shape through to biases on cosmological parameters. In Sec- 
tion [3] we describe the separate contributions to the PSF 
model and in Section |4] we predict the range in PSF sizes 
expected for two different broad-band filters using a realis- 
tic distribution of stellar and galaxy SEDs and investigate 
two different methods for reducing the bias. These involve 
using: (i) galaxy and stellar colours to correct the PSF size 
and (ii) template fitting. Also we place requirements on a 
simple model for the wavelength dependence for a Euclid- 
like survey. Finally, in Section [S] we discuss our findings. 



2 IMPLICATIONS OF USING AN 
INCORRECT PSF 

We first make a simple calculation of the shear bias induced 
by mis-estimating the size and ellipticity of the PSF and 
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Figure 1. A sample of stellar and galaxy spectra showing the 
variety of SEDs among astronomical objects. An elliptical galaxy 
at 2 = 0.9 (solid red), an irregular galaxy at z = 0.7 (dot-dashed 
blue), an AOV star (dashed green) and a KOIII star (dotted ma- 
genta) . 



compare it to the g enera l systematics limit calculated in 
lAmara fc RefregieJ (|2007l ). We then consider the redshift 
dependence of the shear biases and propagate this through 
to biases on cosmological parameters. 

The shear measurement bias caused by an incorrect 
PSF will depend on the shear measurement method em- 
ployed. For simplicity we assume here that shears are cal- 
culated using unweighted quadrupole moments. Although 
this method is not feasible in practice due to the low signal- 
to-noise level in real images, it is relate d to the widely use d 
Kaiser, Squires and Broadhurst method (jKaiser et al.lll995l l. 
It has the advantage of being extremely easy to use for shear 
measurement bias calculations. 

Shear mis-estimates are often quantified by Taylor ex- 
panding the estimated two-component shear 7i (i = 1, 2) in 
terms of the true shear 7^ as 



ii = rUi'yi + a 



(1) 



where rrii is r eferred to as the mult iplicative bias, d as the 
additive bias mans et al.ll2006h and it is often assumed 
mi ~ m2 = m and ci ~ C2 = c . 

lAmara fc Refregieij (l2007f) showed that, for a full-sky 
survey (2 x 10* square degrees of extragalactic sky) with 
35 galaxies per square arcminute and a median redshift of 
0.9, the shear multiplicative error m must be < 1 x 10"'^ 
to keep systematic biases on cosmological parameters below 
random uncertainties, for a range of possible redshift evo- 
lution scenarios for m. Using the equations in Appendix 1X1 
this translates to a requirement on the PSF size mis-estimate 
SFpsF of 
SFpsF 



-PpSF 

lAmara" 



: m s; 1 X 10" 



(2) 



Refregied (|2007l ) also placed a requirement on the 
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mean square error a^y^ < 10^^ for the same survey. Inter- 
preting this as a requirement on the square of the shear 
measurement additive error c^, and neglecting the subdom- 
inant term in Eq. IA8I this places a requirement on the PSF 
ellipticity mis-estimate (5epsFi of 

5ePSF« ~ 4c ^ 1 X 10"^ (3) 

Because the observed SED of a galaxy changes with red- 
shift, the appropriate PSF will also depend on galaxy red- 
shift. Therefore the PSF biases SFpsF and (Jepsp depend on 
redshift and so do the multiplicative and additive biases m 
and c. If, for example, the impact of the multiplicative and 
additive biases as a function of redshift mimics a particu- 
lar cosmological parameter the requirements may be more 
stringent than the above more approximate calculation. We 
will therefore calculate PSF biases as a function of galaxy 
redshift and insert them into the more detailed calculation 
described below, which propagates the effect into biases on 
cosmological parameters. 

Use of the wrong PSF model will cause the measured 
cosmic shear cross power spectra between redshift bins i and 
j, Cij{£), to differ from the true cosmic shear power spectra 
Cijil). If a particular systematic on the cosmic shear power 
spectrum ACfj(^) = Cij{l) - Ctj{t) is ignored then the bias 
on cosmological parameters 5pa is given by 

where i, j, k, I and /? are summed over. Gov ^C^j{£), Cfc;(^)] 
is the two dimensional covariance matrix between the cross- 
spectra and J- is the Fisher matr i x between the cosmologic al 
parameters (|Huterer et al.|[2006l : lAmara fc Refregieij|200'i1 ). 

In the presence of a redshift dependent multiplicative 
bias, the measured lensing power spectrum can be given in 
terms of the true lensing power spectrum by 

(?5W = Cr,W(l+ m^) (5) 

where m' is the multi plicative bias for red shift bin i, aver- 
aged over all galaxies (jHuterer et al. I l2006l . Eq. 16). 

The impact of additive errors depends to first order 
on the spatial variation of the additive errors. For the 
case of a wavelength dependent PSF this will induce power 
on the scale of the separation between stars of a typical 
colour, which may be pr opagated into cosmology (see also 
iGuzik fc Bernsteinll2005l ). Here we focus on PSF size mis- 
estimates and therefore do not consider additive errors fur- 
ther. 

In this paper we compare the PSF sizes for stars with 
those for galaxies, without considering the galaxy morphol- 
ogy or profile. The equations derived above and in the 
Appendix make it possible to draw significant conclusions 
about the cosmology biases independent of considerations 
about the galaxy light distribution, if all parts of the galaxy 
have the same colour. Our main metric is the difference be- 
tween the FWHM of the stellar and galaxy PSFs. We av- 
erage this over populations of galaxies for various different 
PSF correction schemes. 
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Figure 2. Contribution from each compomcnt of the PSF to the 
image quality (FWHM) as a function of wavelength. The dotted 
(blue) line represents the diffraction component, the long-dashed 
(red) line the CCD modulation transfer function component, the 
short-dashed (green) the achromatic component and the thick 
solid line (magenta) shows the overall result taking into account 
all the previously described components. The solid (black) line 
shows the relative instrumental passband using a wide optical 
filter Fl such as that proposed for Euclid. The dashed (black) 
vertical lines show an alternative, narrower, lensing filter (-F4). 
The background shaded areas show approximately the wavelength 
coverage of the grizy filters. 



3 THE PSF MODEL 

We consider a simple instrument model in which the PSF is 
made up of three circular components, each with a different 
wavelength-dependence. This is reasonable for a space-based 
instrument composed mainly of reflective surfaces, such as 
Euclid. For such an instrument the PSF ellipticity is rel- 
atively insensitive to wavelength and the three main PSF 
components are (i) nearly diffraction limited telescope optics 
giving rise to an Airy disk with size inversely proportional to 
the wavelength (ii) the GGD modulation transfer function 
(MTF) which tends to spread out higher energy photons 
more than lower energy photons and (iii) a wavelength in- 
dependent part such as telescope jitter. We assume for sim- 
plicity that each component is Gaussian with a wavelength- 
dependent size. 

We describe the total size of the PSF by its FWHM, F, 
which is given by the quadratic sum of the FWHM values 
of the three components 

F|sf(A) = F^(A) + F^itf(A) + (6) 

Note that the addition in quadrature works reasonably well 
even if the diffraction limited component is not a Gaussian: 
we find adding FWHMs in quadrature works to better than 
5 per cent accuracy when an Airy disk is convolved with a 
Gaussian of the same FWHM, and improves to better than 
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2 per cent accuracy if the ratio of FWHMs is changed by a 
factor of 4 either way. 

The size of the diffraction limited image is given by 



Fd(A) = 0.154" 



D 



A 



(7) 



,1.2my V7350A/ 

where D is the diameter of the primary mirror. We take the 
contribution from the CCD MTF to be that measured em- 
pirically for an e2v CCD 231-84 (M. Cropper, priv. comm.), 
which is given approximately by 

0.11" - 0.027" (^^] if A 5^ 7OOOA 
Fmtf{X)={ „ (8) 

0.09" - 0.007" if A > 7000^1 

V 7000A / 

Finally, we take the contribution to the PSF size from the 
achromatic component to be 



FWHMj = 0.1 



(9) 



as appropriate for a Euclid-like instrument. 

We plot the three contributions to the PSF image size 
in Fig. [2] for a 1.2 meter primary mirror. The image size 
is dominated by the diffraction limit of the instrument. As- 
suming the PSF contributions from different wavelengths all 
have the same centroid we can calculate the FWHM of the 
composite PSF from the FWHM of each component and 
the transmitted flux S'(A)r(A), where S{X) is the spectral 
energy distribution (SED) of the object and r(A) is the in- 
strumental plus filter response, by 



-FpSF 



/ SjX) T(A) F|sf(A) dA 
/ S{\) r(A) d\ 



(10) 



In this paper we assume T(A) is the instrumental response 
for the case where a wide optical Fl filter (5500-9200A; 
see Fig. [5)) is used for measuring the shapes of galaxies, as 
proposed for the Euclid satellite. 



4 PSF AND COSMOLOGY BIASES 

We use realistic galaxy and stellar populations to quan- 
tify the amount by which the PSF FWHM is likely to 
be mis-estimated. Galaxy SEDs are generated using mock 
catalogues designed to simulate the distribution of red- 
shifts, colours and magnitudes of g alaxies in GOODS-N 
ICowie et alJ |2004':'Wi rth et al.ll2004l ). Template spec t ra are 
taken from IColcman ct al.1 (|l980i rand iKinnev et all (119961 ) 
and intermediate types obtained by li near interpolation o f 
these templates (for further details see lAbdalla et aLll2008h . 
In this paper we include galaxies up to a redshift of 2. Galax- 
ies more distant than this will have small apparent sizes 
and are thus unlikely to be used to measure cosmic shear. 
For instance, for the Cosmos survey, much less than 10% of 
the l ensing usable galaxies have z > 2.0 (|Leauthaud et af] 
[2003). The PSF sizes for stars are estimated using stel- 
lar SEDs from the Bruzual-Persson-Gunn-Stryker (BPGS) 
Spectrophotometric Atlafl The catalogue contains 175 dif- 
ferent SEDs covering a broad range of spectral types. The 
true PSF size for each of the star and galaxy types in the 



|http: / / www.stscl.edu/hst / observatory /cdbs / 1 
cal_cataiogs.html#bruzual-persson-gunn-stryker. 
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Figure 3. Histogram of the predicted PSF sizes (FWHM) for the 
mock galaxy catalogue. The average value of the distribution is 
marked on the plot with an arrow, as well as the sizes for typical 
disk and halo stars G5IV and KOHL 



mock catalogues is then estimated by inserting its SED into 
Eq.Iini 

A histogram of the FWHM distribution for the galaxy 
population is shown in Fig. [3] The PSF size ranges from ap- 
proximately 0.175" to 0.220", has an average of 0.1922" and 
a dispersion of 0.006", or 3% of the mean PSF size. As an 
example, the FWHM of a G5 sub-giant star (taken here as a 
typical disk star) is typically 0.1884 arcsec, whereas a KO gi- 
ant star (typical halo/bulge star) is typically 0.1897 arcsec. 
In a conventional analysis which ignores the wavelength- 
dependence of the PSF, galaxies with small angular sepa- 
rations from the above two example stars will often have 
their shears underestimated. In fact, out of all the stars in 
the BPGS catalogue only a quarter of them (cold K and M 
stars) have PSF estimated-sizes larger than the average of 
the galaxies. The multiplicative shear mis-estimates for the 
example G5 and KO stars are of the order of m ~ 3 x 10~^ 
relative to the average galaxy, in clear disagreement with our 
requirements on the PSF size error (Eq. [T| and therefore we 
need to explore methods for mitigating the effect. Below we 
describe two methods for correcting the PSF wavelength- 
dependence using colour information. 

Fluxes are obtained for each object in the mock cat- 
alogues for the filters -Fl, Y, J and H up to the limiting 
magnitudes 26.25, 24.0, 24.0 and 24.0 respectively (AB mag- 
nitudes, 5cr detections). In addition, we consider different 
scenarios for the complementary ground based photometry. 
We assume here observations in the filters g, r, i, z and y 
with three different depths each, shallow, medium and deep 
(see Table [ij . The fiducial optical depth used in this paper 
(medium) is chosen to correspond to a DES or Pan-STARRS 
type survey with two dedicated telescopes (PS2). Shallow 
and deep correspond to a Pan-STARRS type survey with 
one (PSl) or four (PS4) dedicated telecopes. 
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Figure 4. PSF FWHM versus colour for a realistic distribution of galaxies (black dots) and stars (magenta stars) of several different 
spectral types. The continuous line is a second order polynomial fit to the stellar sequence. The insets show the residuals between the 
galaxy and stellar polynomial FWHM values at a given colour, in units of 10"'^ arcsec. Results are shown using the default parameters 
for the mirror (1.2m) and survey depth (medium) for the colours _F1 — Y (top left), r — Fl (top right) and r — FA (bottom right; F4 as 
the shape measuring filter). The effect of using a larger mirror (1.5m) for r — Fl is also shown (bottom left). Fl and F4 correspond to 
the wavelength ranges 5500-9200A and 6120-8580A respectively. 



Table 1. Pan-STARRS (PSl, PS2 and PS4) and DBS optical 
photometry depths. The values quoted here correspond to AB 
magnitudes of 5a detections. 



Band 


PSl 


PS2 


PS4 


DBS 
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24.66 


25.53 


26.10 


25.35 


r 


24.11 


24.96 


25.80 


24.85 


i 


24.00 


24.80 


25.60 


25.05 


z 


22.98 


23.54 


24.10 


24.65 


y 


21.52 


22.01 


22.50 


22.15 



4.1 Broad-band colour method 

Here we investigate the possibility of estimating the galaxy 
PSF from stars with the same broad-band colour. For this 
to work well two conditions must be met (i) there must be a 
good correlation between a given broad-band colour and the 
PSF size and {ii) this relation must be the same for both 
stars and galaxies, despite their potentially differing spectra. 
To assess the extent to which this applies for our fiducial 
instrument we have plotted the PSF sizes for the galaxy 
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and stellar populations against two different colours: F1~Y, 
which can be fully determined by an instrument containing 
a single optical filter plus the Y infra-red band, such as the 
proposed Euclid design; and r — Fl, which requires r-band 
observations which could come from the ground. The results 
are shown in the top two panels of Fig. |4] 

It is clear from Fig. |4] that there is a positive corre- 
lation between the PSF FWHM and the colour i.e. redder 
objects tend to be larger than bluer ones. This is unsur- 
prising given that the contribution to the PSF size from the 
diffraction limit of the telescope is the dominant component. 
It is also clear that the correlation between the PSF size and 
the r — Fl colour is tighter than with the Fl — Y colour. We 
expect this to occur because the r — Fl colour constrains 
the slope of the SED within the lensing measuring filter Fl, 
whereas the Fl — Y colour constrains the SED slope at red- 
der wavelengths, which are less relevant. This interpretation 
is strengthened by the fact that other colours such as r — z, 
g — z and g — Y , which are also able to constrain the slope 
of the SED over the shape measuring filter, produce corre- 
lations almost as tight as those for r — Fl. 

For the distribution of stellar SEDs used in this pa- 
per, we can see from Fig. |4] that the stars occupy a well 
defined locus in the FWHM versus colour space. We fit a 
parabola to this data and by doing so we can estimate the 
PSF size from a colour alone. From the figure we see that 
this line fits well through the stellar data points and thus 
we meet condition (i) for the stars. It is fortunate that the 
wavelength-dependence of the PSF tends to be, to a first ap- 
proximation, relatively stable. This means that it is viable to 
use data from several high signal-to-noise stars from several 
different fields to empirically find this locus (provided that 
the wavelength independent effects on the PSF are properly 
dealt with). To test condition {ii) we have to assess how 
well the PSF size-wavelength-dependence for galaxies fol- 
lows that of the stars. We therefore calculate the residuals 
between the FWHM of the galaxies and the fit to the stellar 
FWHM-colour relation (shown in the insets of the figure) . 

The bias on the PSF size, ((J-fpsp), is reduced dra- 
matically by including ground-based photometry, decreas- 
ing in magnitude from 3.6 x 10~^ arcsec for Fl — Y to 
0.16 X 10^"^ arcsec for r — Fl 0. For a simple calculation 
in which the dependence of the multiplicative bias on red- 
shift is not taken into account, we find from Eq. [5] that 
SFpsY must be < 0.2 x 10~^ arcsec for a typical PSF size 
oi F = 0.2 arcsec. The inclusion of ground-based photome- 
try thus allows us to meet our simple 'back-of-the-envelope' 
requirement (Eq. [2]). Using a shallower r-band photometry 
depth (PSl-like) increases the bias by a factor of 2.5 from 
the value with the fiducial (medium) depth. However, using 
a deeper ground-based photometry depth (PS4-like) has no 
effect on the bias. This suggests that there is an intrinsic 
difference between the distribution of PSF sizes for typical 
star and galaxy SEDs (measured in the Fl band) with the 
same r — Fl colour, which is not reduced by increasing the 
photometry depth beyond 'medium'. 

The wavelength-dependence of the PSF size can be re- 



° Unless otherwise stated, the photometry depths are for a 
'medium-depth' survey, corresponding to a PS2 survey (see Ta- 
ble[T]for depths in the g, r, i, z and y bands). 



duced by increasing the size of the primary mirror (since 
the telescope is diffraction- limited) . We find that the bias is 
reduced by a factor of 1.6 by increasing the mirror diame- 
ter from the fiducial value (1.2m) to 1.5m (Fig. |4] bottom- 
left paneljfl. For instance, for a space telescope, such as 
the HST, with a 2.4m mirror, the diffraction component 
will contribute the same as the CCD MTF and the jitter 
at ~ 7500 A. In this case the overall optics becomes much 
less chromatic, in particular for bluer wavelengths where the 
CCD MTF partially compensates the diffraction effect. This 
is directly reflected in the bias we estimate as (5-Fpsf) gets 
to be as small as 0.03 x 10"'^ arcsec (less than 5 times the 
value of our default configuration). On the other hand a 
mirror as small as 0.6m will induce a bias in the PSF size 
of 0.46 X 10^^ arcsec, or 2.5 times larger than the bias we 
obtained with the fiducial configuration. 

Another way to reduce the chromaticity of the system is 
by using a shape measurment filter F4 whose passband is 2/3 
narrower than Fl (See Fig. [2} .By using the narrower shape 
measurement filter we increased to absoute bias value from 
-f0.16 to -0.18 miliarcsec. As we can see in Fig. |3] (bottom- 
right panel), this configuration indeed decreases the slope 
of the PSF FWHM-colour relation for stars and galaxies. 
However, for this particular colour of choice (r — FA), the 
'matching-up' of stars and galaxies becomes poorer. If, for 
instance, we use the colour g — z instead we get a bias of 
-0.11, that captures the improvement we expected by using 
a narrower filter. 



4-1.1 Redshift dependence 

Galaxies with the same intrinsic spectra will have different 
observed colours as a result of the range of galaxy redshifts. 
The bias on the PSF FWHM is thus redshift-dependent, and 
could potentially be very damaging to weak-lensing tomog- 
raphy. In the previous section we calculated the mean bias 
on the PSF FWHM for the galaxy SEDs averaged over all 
redshifts included in the catalogue [z < 2). In this section 
we investigate the effect of the redshift dependence of the 
galaxy colours. 

Fig [5] shows the redshift dependence of 5FpsF for the 
fiducial scenario (dashed red line). There is a large negative 
bias at high redshift which will affect a small fraction of 
galaxies. We also see clearly the biggest limtation of our 
global average method, since in that method positive and 
negative contributions will cancel. 

We propagate this through to dark energy related cos- 
mological parameters as explained in Section [2] and divide 
the biases on the parameter values by the statistical errors 
on the parameters found using the standard Fisher matrix 
approach, for a 20,000 square degree survey with 35 galaxies 
per square arcminute and a total uncertainty on each shear 
component of a-y = 0.35. The results are shown in Table [2] 

We see that, for all configurations all the dark energy 
parameter biases we get are smaller that the expected sta- 
tistical errors. In particular, for Wa, the most demanding 
parameter those surveys are trying to determine, the bias 

^ For these simulations the changes in the instrumental configu- 
ration do not change the signal-to-noisc of the observations, thus 
the same galaxies are observed in all cases. 
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Table 2. Biases on dark energy related cosmological parameters divided by the statistical uncertainty (0.035, 0.045 and 0.149 for Hde, 
wo and Wa respectively) on the cosmological parameter. Results are shown for each correction method. 
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Figure 5. Residual between the actual and the best estimate 
galaxy PSF FWHM as a function of redshift with the broad-band 
method using the r — Fl colour (red dashed) and the template- 
fitting method using the shape measurement filters Fl(5500— 
9200A) (black solid) and F4 (6120-8580A) (blue dot-dashed). 

corrected by the broad-band colour method is less than a 
half of the statistical error level for our fiducial survey con- 
figuration. The other trends follow those already discussed 
in the context of the global averaged biases discussed in the 
beginning of this Section. 

4.2 Template-fitting inetiiod 

The previous section assumes that we use a single colour to 
determine the correct PSF model to use for a given galaxy. 
In practice we will have more than two filters, which will 
be used to calculate photometric redshifts. We therefore 
also consider the use of a template fitting method to pre- 
dict the PSF FWHM of a galaxy by using all the colours 
available. Using ANNz (jCollister fc Lahavlliooi ) we trained 



and validated a neural network using approximately one 
third of the simulated gala:xies available, to predict the red- 
shifts, spectral types and reddening of each galaxy, given 
the multi-colour information. With this information we can 
compute the SED of each object and use a model of the PSF 
wavelength-dependence to predict the PSF FWHM for this 
galaxy. 

A telescope model and stars will be used to build this 
model for the wavelength-dependence, and the accuracy of 
the model will depend on the stability of the wavelength- 
dependence with telescope properties, and the number of 
stars available to calibrate which model to use. In the case 
of the Hubble Space Telescope, a PSF model taken from 
the telescope design is routinely used in conjunction with 
calibration from any stars in the field to assess the telescope 
configuration in a given observation. Therefore in this paper 
we use the exact model as given in Eq.[T2]and propagate the 
noisy and potentially biased galaxy SED estimates through 
to PSF biases and cosmological parameter biases. 

The comparison between this predicted PSF FWHM 
and the truth for the exact galaxy SED and redshift can be 
seen in Figure [6] We can also compare the global averaged 
bias to compare with the results we got from the broad-band 
method and understadn the trends. By using our standard 
configuration we obtain a bias of {(5Fpsf) = +0.19 miliarc- 
sec, which is similar from what we got by using one single 
colour and also met our back-of-envelope requirment (See 
Fig. E]). The use of deeper photometry, a larger 1.5m mirror 
or the narrower lensing filter F4 reduced (JFpsf) to -1-0.18, 
+0.13 and +0.12 respectively, following the expected trends. 

The redshift dependence of the PSF FWHM bias is 
shown for the fiducial scenario (solid black line) and for a 
senario with the lensing filter F4 (dot-dashed blue line) in 
Fig. [5] Both cases show much less redshift evolution than 
the broad-band colour method, in particular for z > 1.0. 

We propagate this (and all the other scenarios) through 
into biases on cosmological parameters and find the results 
given in the second line of Table [2] All the biases are smaller 
than the statistical errors. In Fig. [7] we show the statistical 
confidence level uncertainties in the space wq — Wa for our 
fiducial scenario. It can be seen in the figure that the resid- 
ual bias after the correction of the wavelength-dependence 
is well within the 68% CL. In comparison to the broad- 
band colour method the template fitting method produces 
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Figure 6. Comparison between tlie actual PSF sizes for objects 
with the SED of the galaxies of the mock catalogue and calculated 
PSF sizes given redshifts, reddening and spectral types obtained 
through ANNz. 
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Figure 7. Bias in the dark energy equation of state parameters. 
The eUipse shows the region of the wq — Wa space contained in 
the 68% confidence level contour for a weak-lensing Euclid-type 
survey, where the central values are respectively 0.0 and -0.95 
(marked by a cross). The asterisk shows the effect of the residual 
bias induced by the wavelength-dependence effect. 



smaller biases on cosmological parameters and takes more 
advantage of the deeper photometry. In this case by using 
the deepest photometry (PS4-like) we could see actual im- 
provement when compared to the default, medium depth 
(PS2/DES-Uke). 



4.3 Requirements on a simple 

wavelength-dependence model 

Now we consider requirements on the parameters of a 
simple model for PSF FWHM wavelength-dependence for 
the template fitting method. We approximate the FWHM- 
wavelength relation to be a simple linear function. The fidu- 
cial configuration is most closely approximated by a straight 
line whose slope is equal to 1.63 x 10~^ arcsec/A with 
FWHM(7350A) = 0.192". We consider a range of slopes, 
from a pure wavelength independent PSF to a relation twice 
as steep as the fiducial configuration. In all the cases we kept 
the PSF FWHM constant at A = 7350A, which are the cen- 
tral wavelenth of the lensing filters we are considering. In 
Figure |5] we show the variation of the systematic to statis- 
tical error for the dark energy parameters Qde, wo and Wa 
for the case where the wavelength-dependence of the PSF 
has been corrected using the template-fitting method. 

In this Figure, as expected, one can see that the bias 
on cosmological parameters increases as the wavelength- 
dependence of the PSF gets stronger. For all dak energy pa- 
rameters the biases lie within the region where the statistical 
errors are larger than the systematic ones \b{p) /a{p)\ < 1. 
For the evolution of the dark energy equation of state pa- 
rameter the normalised bias \b{'Wa)/cr{wa)\ approaches unity 
when the slope is double the value for the fiducial telescope 
mirror size and filter width. 
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Figure 8. Variation of the systematic to statistical error ratio 
(b{p)/<j(p)) of dark energy parameters as a function of the slope 
of the PSF FWHM wavelength relation. The dashed (green), dot- 
dashed (blue) and solid (red) lines represent the error ratio for 
the parameters CdEi '^'0 and Wa, respectively. 
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5 DISCUSSION 

The wavelength-dependence of the point spread function is 
an effect that has to be carefuUy considered for the next 
generation of cosmic shear experiments, particularly if wide 
bands are used for imaging. The different spectral energy 
distributions of stars and galaxies means that the point 
spread function obtained from stars is not the same as that 
for the galaxies and this can lead to a non-negligible bias in 
shear measurements. 

We have, for the first time, set out a formalism for test- 
ing the wavelength-dependence of the point spread function 
for diffraction limited imaging, using the parameters of a 
Euclid-like survey. Given these characteristics the dominant 
wavelength-dependence comes from the PSF size, which we 
parameterise by the PSF full-width at half maximum inten- 
sity. We have shown that the fractional difference in PSF 
FWHM between the stars and the galaxies must be smaller 
than 1 X 10"^. We find the same fractional requirement on 
the PSF ellipticity difference. We have illustrated the for- 
malism on a fiducial Euclid-like telescope for which only the 
PSF size is significantly wavelength dependent. 

We investigated two different methods for correcting the 
effect and found that although they give similar results for 
the average PSF error, the different dependencies of the PSF 
size bias on redshift leads to very different implications for 
biases on cosmological parameters. For this type of analysis 
it is therefore necessary to take into account the redshift- 
ing of galaxy spectra and the cosmological parameters of 
interest. 

The first correction method we consider matches stars 
of a given colour to galaxies of the same colour. This is 
not expected to be a perfect correction method because 
the spectral energy distributions of two objects with the 
same colour are difi'erent. The stars and galaxies need to be 
matched up so they have similar SEDs within the imaging 
band used for cosmic shear galaxy shape measurement. This 
is best achieved if the colour considered matches well to the 
imaging band. We find that a telescope with a wide optical 
band plus infra-red bands cannot sufficiently self-correct for 
PSF wavelength-dependence using the optical minus infra- 
red colour {Fl — Y). This is expected because the difference 
in luminosity between bands so widely spaced is not well 
correlated with the variation in luminosity within the opti- 
cal band itself. The PSF FWHM of the stars is on average 
3.6 X 10~ arcsec and a fractional error of 18 X 10"^, much 
greater than our requirement (1 x 10"'^) . 

To realise the full potential of cosmic shear, all planned 
surveys will estimate the galaxy redshift using photometric 
redshifts. This places stringent requirements on having ad- 
ditional photometry in multiple wavebands, which may be 
obtained from the ground or from space. This provides the 
ideal input into correction for the wavelength-dependence of 
the PSF. We consider the colour r— Fl, which provides much 
more useful information about the spectral energy distribu- 
tion in the Fl optical band. The average difference in PSF 
FWHM between the stars and galaxies now meets our re- 
quirement. The redshift dependence shows some evolution 
and tends to present larger values for redshifts larger than 
1.0. However we find that the dark energy equation of state 
evolution parameter will not be biased by more than the sta- 



tistical error if this method is used, even in the case where 
we use a shallow ground based photometry. 

We therefore consider an additional method in which 
the full range of available wavebands are used, to match 
those used in photometric redshift estimation. We take ad- 
vantage of the fact that many methods for estimating the 
galaxy photometric redshift also provide an estimate of the 
galaxy spectral energy distribution. This allows a full model 
of the galaxy spectrum to be used in correcting for the PSF 
wavelength-dependence. The extent to which this is help- 
ful depends on how well the instrument PSF wavelength- 
dependence is already known. We do already have a reason- 
able model for the wavelength-dependence. Additionally we 
expect the wavelength dependence to be quite stable as a 
function of time, and furthermore it should be very well cal- 
ibratable using stars. In this work we assume that the model 
for the wavelength-dependence is therefore very well known, 
and the limiting factor in the analysis comes from uncertain- 
ties in the galaxy photometry, which limit our knowledge of 
the galaxy redshift and spectral energy distribution. We find 
for the several cases we tested that the redshift dependence 
of the bias is smaller, when compared to the broad-band 
method and, as a consequence, the dark energy cosmologi- 
cal parameter biases are also smaller. 

We consider the effect of a different telescope mirror size 
and imaging filter width. We find that the large mirror size 
does decrease the biases, as expected due to the smaller con- 
tribution to the wavelength-dependence from the diffraction 
limit. The reduction is around 20 per cent for both correc- 
tion methods. Using a narrower filter reduces the scatter 
on the PSF FWHM but does not decrease the bias, for our 
particular configuration. 

Finally we consider a general linear PSF wavelength- 
dependence which matches well to the Euclid-like fidu- 
cial configuration. This allows a requirement on the linear 
slope of the PSF to be obtained for a given required ac- 
curacy on cosmological parameters. We find that, by us- 
ing the template fitting method, a survey in which the 
wavelength-dependence of the PSF is twice as strong as our 
fiducial Euclid-like survey just meets the requirement that 
\b{p)/aip)\ < 1. 

In this work we have studied the first order effects 
of a wavelength dependent PSF and have shown that this 
effect can be mitigated with the addition of photometric 
data. The next step is to consider higher order effects such 
as colour gradients and the spatial correlation function of 
galaxy colours. These issues will be tackled and discussed in 
later work. 
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weighted quadrupole moments Qij as 



£2 = 



Qii + Q22 
Qii — Q22 
Qii + Q22 
2Qi2 



(Al) 
(A2) 

^ 4-n ■ (A3) 

It is shown in IPaulin-Henriksson et al] (|2008h that the sys- 
tematic bias on a galaxy ellipticity component ^e^^^j^ can be 
approximated in terms of the mis-estimates of the PSF size 
SRpsF and PSF ellipticity Jepspi as 



/ -Rpsf 



2 (egaii - epSFi) ^^^^^ - Se-psFi ) (A4) 

rtPSF 



where egai ; are the original (pre-PSF but post shear) galaxy 
ellipticity components and _Rgai is the galaxy size. Similar 
definitions apply for the PSF ellipticity and size. This as- 
sumes the PSF and galaxy size and ellipticity measurements 
are made using unweighted quadrupole moments. This prop- 
agates into a bias on the shear estimate 7i = 7i -I- as 



P7 



(A5) 



where is the shear responsivity given by P"' = 2 — (|ej^) ~ 
1.8 for simple shear measurement methods. Similarly 7^ = 

In this paper we quantify object sizes in terms of the 
slightly more intuitive quantity, the Full- Width at Half- 
Maximum (FWHM), F where 



F 



(A6) 



for a Gaussian profile, which we use for the PSF component 
in this paper. We combine Eq. [1] Eq. IA6I and Eq. IA5I and 
compare with Eq. IA4l to find 
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where the second term in q dominates for typical future sur- 
veys because epspi is small, in addition to the small value 
of the fractional uncertainty in the PSF FWHM. We have 
ignored the contribution to d from intrinsic galaxy elliptic- 
ities, as appropriate for the case where the average c; over 
randomly oriented galaxy ellipticities is required. Where use- 
ful we have converted PSF sizes into FWHM values but we 
have left galaxy sizes written using R since the conversion 
from _R to F is dependent o n object profile, and is far from 
Gaussian for galaxies. As in iPaulin-Henriksson et al] (|2008l ) 
we assume -Rgai > 1.57?psf and thus Rpsp/Rgai — 0.5. Thus, 
as expected, multiplicative errors arise from errors in the 
PSF size and additive errors arise from errors in the PSF 
ellipticity. 



APPENDIX A: MULTIPLICATIVE AND 
ADDITIVE SHEAR MEASUREMENT ERRORS 
FROM PSF MIS-ESTIMATES 



We follow IPaulin-Henriksson et alj (|2008l ). who define ob- 
ject size R and two-component ellipticity e in terms of un- 



